Models for the magnetic ac susceptibility of granular superferromagnetic CoFe/Al 2 Oj 
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The magnetization and magnetic ac susceptibility, \ = %' ~ *x"> OI superferromagnetic systems 
are studied by numerical simulations. The Cole-Cole plot, x" vs - x'> i s used as a tool for classifying 
magnetic systems by their dynamical behavior. The simulations of the magnetization hysteresis 
and the ac susceptibility are performed with two approaches for a driven domain wall in random 
media. The studies are motivated by recent experimental results on the interacting nanoparticle 
system Co8oFe2o/Ab03 showing superferromagnetic behavior. Its Cole-Cole plot indicates domain 
wall motion dynamics similarly to a disordered ferromagnet, including pinning and sliding motion. 
With our models we can successfully reproduce the features found in the experimental Cole-Cole 
plots. 

PACS numbers: 75.60.-d, 75.75.+a, 75.40.Gb, 75.40.Mg 



I. INTRODUCTION 

The physics of interacting ferromagnetic (FM) 
nanoparticlcs is a vivid topic of modern magnetism re- 
search. This also applies to the study of the reversal dy- 
namics in thin FM films. The first subject, the properties 
of interacting FM nanoparticles, is investigated by many 
groups focusing either on the preparatio n (e.g. , Ref. QQ) 
or the magnetic properties (e.g. Refs. Hll3 l4l5l laFn) . Nu- 
merous theoretical studies were perfomed in order to un- 
derstand the observed phenomena or to explore possible 
new effects (e.g. Refs lslllllllTlllTl . 

While individual single-domain FM nanoparticles ex- 
hibit superparamagnetic (SPM) behavioriiii^i^, inter- 
acting ensembles lead to very different kinds of phe- 
nomena depending on the type and strength of inter- 
actions. Dipolar interactions become relevant since the 
magnetic moment e.g. for particles with diameter 5nm is 
of the order 5000/is, while the inter-particle distances 
are of the order 1 — 10 nm. The simple formula for 
the mean dipolar energy of a particle to one neighbor, 
Ed-d/ks = (/^o/47rfcB)/^ 2 /£> 3 , yields already 16 K for 
D = 10 nm. Considering many neighbors and shorter 
distances it is obvious, that the effects of dipolar inter- 
action can be observed even at temperatures of the or- 
der 100 K. In addition, several other types of interac- 
tions are proposed, i.e., higher order multipole terms of 
dipolaniSiii, tunneling exchanged or even retarded van 
der Waals interactions— Independent from the still open 
question which interactions arc relevant, one can summa- 
rize, that essentially three different kinds of phenomena 
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occur—: 

For large inter-particle distances, and hence a small 
concentration of particles, the (dipolar) interaction is 
only a perturbation to the individual particle behavior 
and no collective behavior is foundi£*22*2i For intermedi- 
ate concentrations a superspin glass (SSG) phase is en- 
countered. In this case the particle moments (supcrspins) 



collectively freeze into a spin glass phase below a criti- 
cal temperature, T g ,422i22i£i24 For even higher concen- 
trations a superferromagnetic (SFM) state is found. It 
is characterized by a ferromagnetic arrangement of the 
moments— i 26 4 27 J 28 4 29 J 30 rpj ic ma g n& t;i c dynamic behavior- 
resembles at the first glance that of the SSG case, but 
actually shows features of domain wall motion similar to 
an impure ferromagnet as will be discussed below— Here 
one should mention that also additional types of collec- 
tive ordering are proposed in literature, e.g., the corre- 
lated superspin glass state (CSSG)^ 2 ^ 3 - and also that the 
effects of surface spin disorder may become significant— 

The second topic, the reversal dynamics in thin fer- 
romagnetic films finds equally large interest. Both 
experimental 3 ^*^ and theoretical 3 ^* 3 ^* 3 ^* 3 ^ investigations 
are performed in order to achieve a better understanding 
of the processes during the hysteresis cycle. The mag- 
netization reversal occurs either by domain wall (DW) 
nucleation and motion or by magnetization rotation— 
The DW motion at constant (dc) fields is characterized 
by three regions depending on the field strength, that 
is creep, depinning, and sliding motion. Creep is the 
thermally activated motion of DWs, where the average 
DW velocity is v(H) ex cxp[-(r p /r)(g/g p )-^]i 41 i 42 i 43 
This behavior is encountered at small applied fields, 
H <C H p , where H p is the critical depinning threshold 
and T p proportional to a characteristic depinning energy, 
Up = ksTp. At zero temperature a dynamic phase tran- 
sition of second order at H = H p is found. The mean 
DW velocity, v, can be interpreted as order parameter 
of the depinning transition, with v(H) oc (H — Hp) 13 : — 
At T > the phase transition is smeared out and the 
v{H) curve is rounded. Beyond the depinning region, 
H 3> Hp, sliding motion sets in and the DW velocity 
becomes linear with applied field, v ~ jH. Here 7 is the 
mobility coefficient—* 4 ^ 

In alternating (ac) (magnetic) fields, H = Hosm(uit), 
additional dynamical effects will arise. The coercive 



field and the loop area become dependent on the fre- 
quency or in other words on the field sweep rate; 36 ' 39 dy- 
namic phase transitions and crossovers occurf22^£^£ the 
ac susceptibility vs. temperature shows similar features 
as spin glass systems^ and a DW velocity hysteresis is 
found^ 9 - Different models are employed, i.e., numerical 
solutions of the coupled differential equations of the DW 
displacement starting from Maxwell's equations^, using 
an interface depinning model for an clastic DW in ran- 
dom media™4iii22iiiS (sometimes referred to as quenched 
Edwards- Wilkinson^ (EW) equation), kinetic simula- 
tions of a DW in the sliding motion regime 3 ^ and cal- 
culations based on Fatuzzo's domain theoryS applied to 
ultrathin magnetic layers^ 

In this paper we will present model investigations 
motivated by recent experiments on the SFM system 
[CosoFe2o(l-4 nm)/ Al 2 03(3 nm)]i being a realization 
of a densely packed ensemble of interacting nanoparti- 
cles. The complex magnetic ac susceptibility, x' ~ *x" \ 
reveals that the magnetic dynamic behavior can be ex- 
plained within the concept of domain wall motion in an 
impure ferromagnet^*^ 3 " That means, the granular sys- 
tem behaves like a thin FM film, only with the differ- 
ence, that the atomic moments arc to be replaced by 
'super- moments' of the individual particles. This con- 
cept implies that the FM nanoparticles remain single- 
domain whereas the ensemble shows collective SFM be- 
havior. This idea is evidenced from the Cole-Cole plot, 
x" vs. x'"^ Hence we will focus on the Cole-Cole pre- 
sentation and compare it to that found experimentally. 



II. AC SUSCEPTIBILITY AND COLE-COLE 
PLOTS 



Magnetic systems exhibiting relaxational phenomena 
can be characterized by the complex ac susceptibility, 
x(w) = x' ~ ix" ■ The time dependent complex ac sus- 
ceptibility is defined as 



M(t) = x(t)H{t) , 



(1) 



with the complex external field H(t) = —iH^e tu}t [H(i) = 
Rc(H(t)) = H sm(ujt)] and the magnetization M. In 
this paper we study the time independent term of the 
Fourier series for x(t), namely: 



x = x' 
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This defines x' an d x" as follows 
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One way of presenting the data is the Cole-Cole or 
Argand representation. The imaginary part is plot- 
ted against the real part of the susceptibility, x" 
vs. -^',5^55 jj- can serve as a fingerprint to distin- 
guish different magnetic systems by their dynamic re- 
sponse. E.g. a monodisperse ensemble of non-interacting 
SPM particles has exactly one relaxation time, r = 
to exp(KV/ k B T) jlili and will display a semicircle with 
the center on the x'-axis. Here K is an effective 
anisotropy constant, V the volume of the particle and 
To corresponds to the microscopic spin-flip time which 
is of order of 10~ 10 s. The Cole-Cole plot can easily be 
derived from an analytic expression for the ac susceptibil- 
ity given in Ref. |S| for a monodisperse SPM ensemble in 
zero-field with a random distribution of anisotropy axis 
directions: 



x H = ^— 
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(7) 
(8) 



where M s is the saturation value of the magnetization. 
Defining a = /.iqM^/3K and a = KV/k B T and eliminat- 
ing to one gets 



X 




(9) 



with T = 2-k/lu = 1//. 



which describes a circle with the radius r = aa/2 and 
center at (a(2 + c)/2; 0) in the Cole-Cole plane. 

In Fig. [Ha) the result is shown for parameters 
/j, M^/3K = 1 and KV/k B T = 1. In the case of a parti- 
cle size distribution (polydispersivity) and hence a distri- 
bution of relaxation times the Cole-Cole semicircle is ex- 
pected to become flattened and/or distorted^ Fig. IHb) 
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FIG. 1: Cole-Cole plots, \" vs - X j ( a ) analytically obtained 
for a non-interacting monodisperse ensemble of SPM particles 
with hqMs/ZK = 1 and KV/k B T = 1 (see text), (b) numeri- 
cal result for a polydisperse ensemble with a log-normal distri- 
bution (circles) and a Maxwell distribution (diamonds) of par- 
ticle volumes with noM^/SK - 1, K/k B T = 1, r = 1, AV = 
0.9 and (I 7 ) = 1, and (c) shows experimentally obtained 
curves on the SSG system [CoFe(0.9 nm)/Al203(3 nm)]io at 
three different temperatures, T = 50, 55 and 60 K (Ref. l57f) . 
The particle sizes follow a Gaussian distribution with (V) = 
11.5 nm 3 and AV = 0.95 nm 3 . The frequency range is indi- 
cated in the figure. 



shows two numerically obtained curves, where a particle 
volume distribution from a log-normal distribution (cir- 
cles) and a Maxwell distribution^ (diamonds) is assumed 
using /j, M*/3K = 1, K/k B T = 1, r = 1, (V) = 1 and a 
relatively broad distribution width AV = 0.9. One finds 
an asymmetric Cole-Cole plot for the case of a log-normal 
distribution. Obviously this is due to the asymmetry of 
the distribution itself. By choosing the more symmetric 
Maxwell distribution the curve becomes symmetric and 
only slightly shifted downward. Extremely high poly- 
dispcrsivity is found in spin glass systems, where the 
distribution of relaxation times is expected to become 
infinitely broad due to collective behavior^ Fig. ^c) 
shows an experimentally obtained Cole-Cole plot on the 
SSG system [Co 80 Fe 2 o(0.9 nm)/ Al 2 3 (3 nm)] 10 at dif- 
ferent temperatures, T = 50, 55 and 60 KS Here the 
particle sizes follow a relatively narrow Gaussian distri- 
bution with (V) = 11.5 nm 3 and AV" = 0.95 nm 3 as ev- 
idenced from a transmission electron microscopy image 
for a simliar sampled 



III. MODELS 

We study the complex ac susceptibility with two dif- 
ferent approaches, where account is taken of the fact, 
that M is controlled by the field-induced sideways mo- 
tion of one DW. In this case it follows that M{t) oc v(t), 
where v(t) is the (mean) DW velocity, beeing a function 
of the external field H(t) and temperature T. Both ap- 
proaches are based on the same underlying model for a 
(i-dimensional elastic DW in a D = (d + l)-dimcnsional 
random environment: 



n 



d d x 



§(Vx^) 2 



H(t)Z + V R (x,Z) 



(10) 



where Z = Z(~x, t) is the d-dimensional displacement pro- 
file of the DW with internal coordinate x, T the stiffness 
of the DW, and Vr the (quenched) random potential. Vr 
can be written in the following way: 



V«[x,Z(x,i)] = - / dZg[x,Z(yL,t)\ 



(11) 



where g[x, Z(x, t)] describes the random force acting on 
the DW with (g) = and (g(x,z)g(x',z')) = <5 d (x - 
x')Ao(z— z'), with Aq(z) = Ao(— z) being a random force 
correlator which is a monotonically decreasing function 
decaying over a finite distance £. 

Since the experimental system is a magnetic film, we 
restrict ourselves to the case D = 2 in the following. The 
dynamics of the system follows from the EW equation of 
motion: 



ldZ(x,t) 

1 



Of 



SH 
"5Z 



7?(x,i), 



(12) 



where 7 is a kinetic coefficient and ?7(x, t) a thermal noise 
term. The DW velocity is given by u(x, t) = Z(x, t). 
Here we are interested in the mean DW velocity v(t) = 
(i>(x, t)) x and mean displacement Z(t) = (Z(x, t)) x , from 
which we can calculate the ac susceptibility as described 
above. Here (. . .) x denotes the average over the internal 
DW coordinate x. 

(i) adiabatic approach. — We use the expression for the 
mean DW velocity in the adiabatic driving regime follow- 
ing from a functional renormalization group (RG) treat- 
ment of 1)12(1 . given in Ref.|43, which interpolates between 
the creep regime and sliding DW motion, 



where x = H/H p , y = T p /T and 
9(1 -x) 



F(x,y) 



1 + {yx-^y/ 6 
+e(x - 1) 



exp [yx M (l - xf 



(13) 



(14) 



1 



1 + {yx-^)P/ e 



Here <d(x) is the step function, T p ~ T£ 2 L?~ d the typical 
pinning energy on the Larkin length scale L p , H p the 
zero temperature depinning field, and /z, /?, and 6 the 
relevant critical exponents^ which depend on the DW 
dimension d. A time discretization, At, is used which is 
chosen to be much smaller than the period of the driving 
field, At = 10~ 5 T. Then Z(t) is calculated for each time 
step by a simple integration of Eq. ^3 i.e., AZ(tj) = 
v(H(ti))At, where v(H(ti)) is calculated for each time 
step from Eq. El Here the values of time t, T, u), and 
/ are choosen to be dimensionless, since no quantitative 
comparison to the experiment is required. Formally this 
can be done by introducing an arbitrary time scale to and 
substituting t — > t/to. Analogously this can be applied to 
all other parameters and observablcs, i.e. field Hq/H p — > 
Hq, temperature T/T p — ► T, velocity v — > v/(^H p ) and 
lenght, L z /Lq —* L z , where Lq is an arbitrary lenght 
scale. 

The magnetization for a finite system is defined here 
as 



M(t) 



2Z(t) 



(15) 



where L z is the extension of the sample in Z-direction 
and < Z < L z . This implies, that -1 < M < +1. In all 
cases the initial condition is Z(0) = L z /2. This approach 
includes the temperature as a parameter, but we restrict 
our investigations here to small values, T = 0.1. 

(ii) non-adiabatic approach. — Since Eq. I|13|) was ob- 
tained for an adiabatically changing field, it can only be 
used as an approximation, if the frequency is sufficiently 
small^ (see also Fig. 2 in this Ref.). In order to include 
the pronounced non-adiabatic effects at higher frequen- 
cies (e.g., hysteresis of the velocity), one has to start with 
the underlying equation of motion 112|> which yields 



ldZ(pc,t) 

1 



dt 



rV 2 Z(x,t)+iJ(t) + <?(x,Z(x,t)), (16) 



where the thermal noise term is neglected which is justi- 
fied since the relaxation times for the DW creep at low 
temperatures are very long (^> uj^ 1 ) and we consider onl; 
finite (not exponentially small) frequencies. In Ref 
this equation is studied in detail in the case of an ac 
driving force in an infinite system and it is shown that 
thermal effects are not essential for not too low frequen- 
cies (uj > lot ~ up(T/T p y z l with the critical exponents 
v, z, and 0). Therefore we can restrict ourselves to the 
zero temperature equation of motion. 

In this approach we investigate both infinite (L z — > oo) 
and finite (L z < oo) systems. In the second case, the 
DW will hit the boundary of the system for low enough 
frequencies, such that the magnetization will saturate 
(—1 < M < 1). Therefore we can derive a critical fre- 
quency u> c above which the system will behave as an infi- 
nite system. The finite frequency of the driving force acts 
as an infrared cutoff for the propagation of the DW which 
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FIG. 2: M vs. H curves from simulations of the adiabatic 
approach (i) with T = 0.1, H = 1.85, L z = 8.0, n = 0.24, 
8 — 0.83, and /3 — 0.66 at different frequencies / = 1.6 • 10~ 7 
(a), 1.6 • 10~ 3 (b), 7 ■ 10~ 3 (c) and 8 • 10~ 2 (d).' Note that all 
quantities are measured in dimensionless units, as described 
in the text. Lines are guides to the eye. 



can move up to a length scale L w = L p (T"/L p 2 /uj) . 
Equating this scale to L z gives the following expression 
for lo c (f c accordingly) 



LU C « lo p (L p /L z ) 2 



(17) 



with the typical pinning frequency ui p = jT/L p . 

For the numerical integration of Eq. (|16fl it is dis- 
cretized in x-direction(s) into N d positions with a lattice 
constant a. Here we also go over to dimensionless units 
with an arbitrary time scale to. These two parameters, 
a and to, are chosen such that to7r/a 2 = 1 and that the 
dimensionless random force tojg are set to values in the 
interval [—1/2, 1/2] at positions with distance £. Between 
these positions g is interpolated linearly which results in 
a Gaussian distribution Aq(z) with variance I. The de- 
pinning field H p is not used as input parameter but can 
be calculated from Eq. I|16|) (at u> = 0) using a bisection 
procedure with constant amplitude. 

For our simulations we choose £ = 0.1, N — 1000, for 
finite systems L z = 8.0, and a time discretization such 
that At <C min(o;~ 1 , 0.1). The results for \ are averaged 
over 100 disorder configurations for each frequency. 



IV. RESULTS AND DISCUSSION 

Fig. |3 shows an example of hysteresis loops from sim- 
ulations within approach (i) with T = 0.1, Hq = 1.85 
and L z = 8.0 at different frequencies / = 1.6 ■ 10 -7 (a), 
1.6 ■ 10~ 3 (b), 7 ■ 10~ 3 (c) and 8 • 10~ 2 (d). Note that all 
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FIG. 3: (a) ac susceptibility, x' an d x" vs - ac frequency, /, 
obtained with model (i) with same parameters as in Fig. [2] 
(b) Same data, but plotted in the Cole-Cole presentation, x" 
vs. x' ■ The solid line represents a least square-fit of the low- 
frequency data to a circle and the arrow shows in direction of 
increasing /. 



quantities are measured in dimcnsionless units, as men- 
tioned above. For the values of the critical exponents, 
we use the results from the RG for d = 2, i.e., fi = 0.24 
(Ref. E|), 6 = 0.83 (Ref. [H, and /? = 0.66 (Ref. E| 
(Note, that the precise values of these exponents do not 
have a significant influence on the behavior under con- 
sideration here, especially the qualitative picture does 
not change if the values are modified slightly.). With 
increasing frequency, the hysteresis loop broadens until 
it becomes elliptically shaped above / = 10~ 2 loosing 
also its inflection symmetry. Similar results are found in 
cxperimentsi^LSi 

The ac susceptibility of such hysteresis cycles can be 
calculated from equations J2} and (gj . In Fig. |3 the ob- 
tained data is shown for the same set of parameters as 
for Fig. [21 In (a) one finds the real and imaginary part 
of the ac susceptibility, \' and x" > as function of the ac 
frequency. The real part shows an order-parameter like 
behavior with a non-zero value below and a vanishing 
value above f c « 10~ 2 . Furthermore, the imaginary part 
has a peak at / = 8 • 10~ 3 « f c . 

In the Cole-Cole plot, Fig. [3] (b), this transition ap- 
pears as a sharp change of the slope and curvature. At 
low frequencies, / < f c one observes a quarter-circle cen- 
tered on the x' ax is- It is possible to fit a circle with the 
center on the x'-axis to the low-frequency data (see solid 
line in (b)). This corresponds well to the experimental 
result 31,53 (Fig.|2} and suggest the existence 1 of one effec- 
tive relaxation time in the system. However, for / > f c 
only a vertical line can be observed. This result differs 
from that found in experiment, where the high-frequency 
part is characterized by a positive slope and positive cur- 
vature. This discrepancy needs a closer inspection here. 

By comparison of the susceptibility data to the cor- 
responding hysteresis loops (Fig. [2J, one sees, that f c 
marks the transition between loops saturating at high 
fields (low-f) and those, which do not saturate (high-f). 
In the second case the domain wall is always in motion 
throughout the entire field cycle. The real part is then 
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FIG. 4: Experimental Cole-Cole plot taken from Ref. |3l] 

showing x" vs - x' obtained on the SFM granular sys- 
tem [CoFe(1.4nm)/Al203(2nm)]io. The susceptibility was 
mesured at ac amplitudes fioHo = 50 (a) and 5//T (b) at 
10 mHz < / < 1 kHz at T = 380 (1), 350 (2), 320 (3) and 
260 K (4). Transition fields are marked by arrows—. 



zero, whereas the imaginary part has a 1// dependence 
[Fig. |3](a)], which follows directly from our result shown 
in Ref. l3lL where the complex susceptibility in the case 
of sliding DW motion is given by x = Xoo[l + l/ii^r)], 
or more generally by x = x'oo + x'L/( iuJT )- For x'oc = ° 
this yields directly the vertical part in the Cole-Cole plot 
[Fig. 01(b)]. In Ref. l3ll was argued, that the non-linearity 
of the v(H) function in the creep regime can be taken into 
account by introducing a polydispersivity exponent,/?, in 
the above equation, x = Xoc[l + l/(* w ' r )' 3 ] (compare 
to a similar relationship formulated for the conductiv- 
ity of disordered hopping conductors <j(w) ~ (— iuJT) v ^ T > , 
where < v < 1)SLS1. This yields the linear relationship, 
x" = tan(7r/3/2)[x' — Xoo]- Note, that for any velocity 
function v = v(H) with v(H) = —v{—H) and without ve- 
locity hysteresis*^ it follows, that x' = an d x" ^ l//^ 1 
This can easily be seen from Eqs. (J5) and JHJ and M oc v. 
The consequence is that a monotonically increasing part 
with finite slope cannot be found in the Cole-Cole plot 
by considering only the adiabatic motion of one DW. 

There are two possible ways to improve the model. The 
first one is to simulate an ensemble of non-interacting 
subsystems with different domain propagations lengths, 
pinning fields, H p , or depinning energies, T p . It is possi- 
ble, that this case would yield the situation above qualita- 
tively described by the polydispersivity exponent (3 < 1. 
The second is to employ a more realistic description of 
the DW by using the above introduced non-adiabatic ap- 
proach (ii). The latter case was studied here. 

In Fig. the results for the magnetization hysteresis 
of a DW from Eq. ^Jfor Hq = 1.85 are presented. The 
plots (a) to (c) show hysteresis loops at different frequen- 
cies, / = 0.0016 (a), 0.08 (b), and 0.48 (c) for an infinite 
system (L z —* oo). Here we define M = Z. In this 
case the DW never touches the sample boundary. At 
low frequencies one finds a symmetric loop with respect 
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FIG. 5: M vs. H curves for the non-adiabatic approach (ii) 
with Ho = 1.85 at different frequencies, / = 0.0016 (a), 0.08 
(b), and 0.48 (c) for an infinite system (in this case one defines 
M — Z). (d) shows the magnetization curve for / = 0.0016 
but a finite system (L z = 8.0) so that the DW touches the 
boundaries. 
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FIG. 6: Real and Imaginary part of the ac susceptibility vs. 
frequency, calculated within approach (ii) for infinite systems 
(a) [the real part is shifted and scaled] and the corresponding 
Cole-Cole plot (b). In (c) \' an d x" are plotted for the finite 
system (L z = 8.0) and Cole-Cole plot (d). The inset in (d) 
shows the high frequency behavior in more detail. The arrows 
show in direction of increasing frequencies. All simulations 
were performed with Ho — 1.85. 



to the M axis (a) similar to the result shown above in 
Fig- Eld). This symmetry is lost upon increasing the fre- 
quency, (b) and (c), and the loop becomes tilted. This 
tilting is responsible for a non- vanishing real part of the 
ac susceptibility and cannot be observed in approach (i). 
The tilting corresponds to the appearance of a velocity 
hysteresis^ That means there exists no functional rela- 
tionship between the velocity and the field any more, as 
it is the case in the adiabatic regime. 

The resulting susceptibilities are plotted in Fig. |5] In 
(a) and (b) the real and imaginary part vs. log(/) and 
the corresponding Cole-Cole plot, respectively, are shown 
for an infinite system, when the DW never touches the 



boundary. In (c) and (d) the same plots are shown for a fi- 
nite system (L z = 8.0). While the low-frequency parts re- 
semble those from approach (i), the high-frequency part 
shows a completely different behavior. For x' -> we 
find in the Cole-Cole plot [inset in Fig. 0(d)] a curve with 
positive curvature similar as in the experiment (Fig.^J. 
One can expect that x goes to with u> — > oo, since the 
velocity hysteresis disappears for u) — > oo. Obviously the 
more realistic second model is capable to describe the ex- 
perimentally found behavior. At this point we want to 
emphasize that the adiabatic approach only works for low 
frequencies, where non-adiabatic effects can be neglected. 
Furthermore it only works at finite temperatures. On the 
other hand the non-adiabatic approach can explain the 
main experimental features even if we use the zero tem- 
perature equation of motion, since the smearing effects 
of the depinning transition due to finite frequencies dom- 
inate the thermal creep effects at low temperatures. 

However, two drawbacks still exist. One, the Cole-Cole 
plot from the simulation shows a rather steep and nar- 
row increasing part compared to the experiment. Second, 
we cannot retrieve the experimentally observed saturat- 
ing part for the highest experimental frequencies, where 
the imaginary part becomes constant (see Fig.0] inset). 
This case was attributed^! to the reversible relaxation 
response of the DW for high frequencies and small exci- 
tations fields i 50 1 58 It would be interesting to study this 
case with a suitably modified model which includes mul- 
tiple and interacting DWs. 



V. CONCLUSION 

In order to get a better understanding of the mag- 
netic behavior found in the superferromagnetic granu- 
lar multilayer [Co8oFe2o(l-4nm)/Al203(3nm)]io we em- 
ployed two types of simulations of a domain wall in ran- 
dom media driven by an external magnetic field. Us- 
ing the first approach with the mean velocity of a do- 
main wall in the adiabatic limit, one can explain the 
monodisperse dynamic response evidenced by a partial 
semicircle centered around the \' axis. However, it fails 
to describe the increasing part with positive curvature for 
higher frequencies in the Cole-Cole plot. This behavior 
can be found by taking the full equation of motion into 
account, where an elastic interface is driven in general 
non-adiabatically in a random medium. Hence a model 
of an impure ferromagnct is capable to describe the main 
features of the experimental results. We find that the ap- 
pearance of a velocity hysteresis is a crucial element in 
the dynamic response of the superferromagnet. We show 
that a Cole-Cole plot may be used to classify magnetic 
systems by their dynamic response. E.g. the above men- 
tioned granular superferromagnet can unambiguously be 
distinguished from a superparamagnet and a supcrspin 
glass system. 
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